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Abstract 

We give an algorithm that constructs the Hasse diagram of the 
face lattice of a convex polytope P from its vertex-facet incidences in 
time 0(min{n, m} ■ a ■ ip), where n is the number of vertices, m is the 
number of facets, a is the number of vertex- facet incidences, and if is 
the total number of faces of P. This improves results of Fukuda and 
Rosta [5], who described an algorithm for enumerating all faces of a 
d-polytope in 0(min{n,m} • d- ip^) steps. For simple or simplicial d- 
polytopes our algorithm can be specialized to run in time 0(d ■ a ■ (p). 
Furthermore, applications of the algorithm to other atomic lattices are 
discussed, e.g., to face lattices of oriented matroids. 

Keywords: polytope, face lattice, enumeration, vertex-facet incidences, 
algorithm, oriented matroid 
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1 Introduction 

Let P be a d-polytope, i.e., a d-dimensional bounded convex polyhedron. It 
is well-known that the set J- of its faces (including and P itself), ordered 
by inclusion, is a graded, atomic, and coatomic lattice: the face lattice of P. 
In particular, each face can be identified with its set of vertices or the set 



*TU Berlin, Fakultat II, Institut fiir Mathematik, MA 6-2, Strafie des 17. Juni 136, 
10623 Berlin, {kaibel, pfetschj@math.tu-berlin.de. The first author has been sup- 
ported by a DFG Gerhard-Hess-Forschungsforderungspreis donated to Giinter M. Ziegler 
(Zi 475/2-3). 



1 



of facets it is contained in. In this paper, a face is usually identified with 

its vertex set. We define ip := \T\ and denote by C the Hasse diagram (as 
an abstract graph) of the face lattice. Hence, £ is a directed rooted acyclic 
graph whose nodes correspond to the elements of If In, are nodes in C 
and H,G e are the corresponding faces of P, then there is an arc {in, (-g) 
in C if and only ii H CI G and dim(G) = dim.[H) + 1. 

The combinatorial face lattice enumeration problem is the following: 
given a vertex- facet incidence matrix of P (see Section 2 for a definition), 
construct the Hasse diagram C of the face lattice. By definition, £ is unla- 
beled. Nevertheless, it might be desired to label each node of C correspond- 
ing to a face F with the set of (indices of) vertices contained in F, the set 
of (indices of) facets containing F, or with the dimension of F. 

Pukuda and Rosta [5] gave an algorithm for the combinatorial face lattice 
enumeration problem for d-polytopcs P which runs in 0(min{n, m} • d ■ ip"^) 
time, where m is the number of facets and n is the number of vertices of P. 
Since (p can be exponential in n and m (consider the d-simplex, for instance) 
it is desirable to have an algorithm whose running time depends only linearly 
on if (and polynomially on n and m). The main result of this paper is such 
an algorithm. 

For the geometric face lattice enumeration problem, which asks for the 
face lattice of a polytope that is given by an inequality description, there are 
algorithms satisfying this requirement on the running time, e.g., by Pukuda, 
Liebling, and Margot [4]. However, in our context no geometric data are 
available. 

Ganter [7] described an algorithm which, given the incidences of atoms 
and coatoms of a general atomic and coatomic lattice, enumerates all elem- 
ents of the lattice in lexicographic order, where each clement is identified 
with the set of atoms below it (which are ordered arbitrarily). Specialized 
to our situation, one obtains an algorithm that computes all vertex sets of 
faces of P in 0(min{n, m} ■ a ■ ip) steps, where a is the number of vertex- 
facet incidences of P. Note that d-max{n, m} < a < n-m,, in particular, a is 
bounded polynomially in n and m. This algorithm, however, does not com- 
pute the inclusion relations between the faces, i.e., the edges of the Hasse 
diagram of the face lattice. Of course, once all (vertex sets of) faces are 
computed, one may construct the Hasse diagram in an obvious way after- 
wards, but this would require a number of steps which is quadratic in the 
total number (p of faces. 

Inspired by Ganter's algorithm, we developed the (quite different) algo- 
rithm presented below, which computes the entire Hasse diagram in the 
same running time of 0(min{n, m} • a • </?), see Theorem 1. It requires 
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0((^ • min{n,m}) memory (without output storage). In our algorithm, the 
vertex set of each face or the set of facets it is contained in, as well as its 
dimension, is readily available (or can be computed without increasing the 
asymptotic running time). Of course, this may increase the (output) storage 
requirements significantly. 

Pukuda and Rosta [5] also considered the combinatorial face lattice enu- 
meration problem for the special case of simple or simplicial polytopes. They 
presented an algorithm that computes the face lattice of a simple polytope 
in 0{d ■ if) steps, provided that in addition to the vertex-facet incidences an 
acyclic orientation of the graph of the polytope is given that induces pre- 
cisely one sink on every non-empty face. Such an orientation is called a good 
orientation or an abstract objective function orientation. Unfortunately, no 
polynomial time algorithm is known that computes a good orientation of a 
simple polytope P — neither if P is given by its vertex-facet incidences nor 
if it is specified by its whole face lattice. 

For simple or simplicial polytopes, our algorithm can be specialized such 
that it computes the Hasse diagram of the face lattice in 0{d ■ a ■ (p) steps 
from the vertex-facet incidences, where no good orientation is required (see 
Section 3.1). 

In Section 2.1 we give a rough sketch of the algorithm, which is followed 
by a more detailed description in Sections 2.2, 2.3, and 2.4. In Section 2.5 
we analyze the algorithm. We present the specialization of the algorithm for 
simple or simplicial polytopes in Section 3.1 and a variant that computes 
the fc-skeleton in Section 3.2. Furthermore, in Section 3.3 a version that 
needs significantly less memory is described which enumerates just the faces 
together with their dimensions (i.e., without the edges of the Hasse diagram). 
Finally, a modification that computes the face lattice of an oriented matroid 
from its cocircuits (Section 3.4) is explained. 

For the basic properties of polytopes that are important in our context, 
we refer to Ziegler's book [8]. The few concepts from the theory of algo- 
rithms and data structures that play a role in the paper can be found in any 
corresponding textbook (e.g. in the one by Gormen, Leiserson, Rivest, and 
Stein [2]). Our running time estimates refer to the uniform time measure 
(i.e., every arithmetic operation/comparison takes one unit of time), while 
our statements on memory requirements refer to the bit model. 
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2 The Algorithm 



Define m to be the number of facets and n the number of vertices of the 
d-polytope P. Let A = {afy) G {0, i}"^^^'^ be a vertex-facet incidence matrix 
of P. Hence the facets of P can be identified with F := {1, . . . ,m} and 
its vertices can be identified with V := {l,...,n}, such that aj^, = 1 if 
facet / contains vertex v, and o/^ = otherwise. Denote by a the number 
of vertex-facet incidences, i.e., the number of ones in A. For S <^V , define 
F(S') := {/ G F : afg = 1 for ah s e 5}, the set of facets containing all 
vertices of S. For T C F, define V(r) := {v e V : atv = I for all t G T}, 
the set of vertices contained in all facets of T. 

For S QV, the set c\{S) := V(F(5)) is the (vertex set of) the smallest 
face of P containing S (in lattice theoretic terms, the join of the elements 
in S). One can check easily that this defines a closure map on the subsets 
of F, i.e., for all 5, S" C 1/ we have: 

S C cl(5), cl(cl(5)) = cl(5), S^S' ^ c\{S) C cl(S'). 

The faces of P correspond exactly to the closed sets of V with respect to 
this closure map (i.e., sets S C.V with cl(5) = S). Our algorithm crucially 
relies on the fact that closures can be computed fast (see Section 2.2). 

2.1 The Skeleton of the Algorithm 

The strategy is to build up the Hasse diagram C of the face lattice from 
bottom (0) to top (P). Consequently, C is initialized with the single face 
and then enlarged iterativcly by adding out-neighbors of nodes that have 
already been constructed. Wc will say that a face has been seen, once its 
corresponding node in C has been constructed. 

During the algorithm, we keep a set Q containing those faces that we 
have seen so far, but for which we have not yet inserted their out-arcs into 
the Hasse diagram. At each major step, we remove a face H from the set Q 
and construct the set Q of all faces G with H C. G and dim(G) = dim(if)-|-l. 
For each face G E Q we check whether it has already been seen. If this is 
not the case, then a new node in C representing G is constructed, and G is 
added to Q. In any case, an arc from the node corresponding to H to the 
node corresponding to G is inserted into C 

In order to compute the set Q, we exploit the fact that Q consists of 
the inclusion minimal faces among the ones that properly contain H. Since 
the face lattice of a polytope is atomic, each face G (z G must be of the 
form H{v) := cl{H U {v}) for some vertex (atom) v; in particular, the Hasse 
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diagram has at most n ■ (p arcs. Thus, we first construct the collection Ti of 
all sets H{v), v eV\H, and then compute Q as the set of inclusion minimal 
sets of H. 

Computing H(v) for some v & V \ H requires determining a closure. 
In Section 2.2, we describe a method to perform this task in 0(a) steps. 
Determining the inclusion minimal sets in the collection Ti clearly could 
be done in 0(n'^) steps by pairwise comparisons, since Ti has at most n 
elements, each of size at most n. In Section 2.3 we show that this can even 
be performed in 0(77,^) time. 

Another crucial ingredient is a data structure, described in Section 2.4, 
that allows us to locate the node in C representing a given face G or to 
assert that G has not yet been seen. This can be performed in 0(a) steps. 

A summary of the analysis of the time complexity of the algorithm, along 
with a pseudo-code description of it, is given in Section 2.5. 

2.2 Computing Closures 

In order to be able to compute closures fast, we store the incidence matrix A 
in a sorted sparse format both in a row and column based way. For each 
vertex v £ V, the elements in F({v}) C {!,... ,m} are stored increasingly 
in a list. Similarly, for each facet / G -F, we store the sorted set V({/}) 
in a list. This preprocessing can be performed in 0{n ■ m) time (which 
is dominated by 0(n • a) and thus does not influence the estimate of the 
asymptotic running time in Proposition 1 below) . The sorted sparse format 
uses 0(a • logmax{n,m}) storage. 

Whenever we want to compute the closure of a set C y, the first step 
is to compute F(5), i.e., the intersection of the lists F({u}), f G 5. Since 
the intersection of two sorted lists can be computed in time proportional to 
the sum of the lengths of the two lists and because the intersection of two 
lists is at most as long as the shorter one, F(S') can be computed in time 
0(^^g_5|F({u})|) C 0(a). Similarly, V(r) can be computed in time 0(a) 
for a set T C F. 

Lemma 1. The closure c\{S) of a set S C.V can he computed in 0{a) steps 
(provided that the vertex-facet incidence matrix is given in the sorted sparse 
format). 

2.3 Identifying the Minimal Sets 

Suppose that H C. V Is a, face of P and Ti. is the collection of all faces 
H{v) = d{H U {v}) <^V,v eV\H. 
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Our procedure to identify the set Q of minimal sets in the collection H 

starts by assigning a label candidate to each vertex in V\H. Subsequently, 
the label candidate of each vertex will either be removed or replaced by a 
label minimal. We keep the following three invariants: For each vertex v 
that is labeled minimal we have H{v) G Q; if two different vertices v and w 
both are labeled minimal, then we have H(v) ^ H(w); Q is contained in the 
set of all H{v) for which v is labeled minimal or candidate. Clearly, if no 
vertex is labeled candidate anymore, the set of vertices labeled minimal is 
in one-to-one correspondence to Q via H{-). 

Suppose there is still some v labeled candidate available. If H{v) \ {w } 
contains some vertex then we have H{w) C H{v), because H{w) is the 
intersection of all faces containing H and w, and one of these faces is H{v). 
Hence, if w is labeled minimal or candidate, we remove the label candidate 
from V] otherwise we label v minimal. 

It follows by induction that the three invariants are satisfied throughout 
the procedure. Moreover, at each major step (choosing a candidate v) the 
number of candidate labels decreases by one. Since each such step takes 
0(n) time, the entire procedure has complexity O(n^). 

Lemma 2. The set Q of inclusion minimal sets inli. = {H{v) : v E:V\ H} 
can be identified in O(n^) steps. 

2.4 Locating Nodes 

During the algorithm, we have to keep track of the faces that we have 
seen so far and their corresponding nodes in C To this end, we main- 
tain a special data structure, the face tree. In this data structure, a face 
S = {si, . . . , Sk} C V (with ,si < • • • < Sk) is represented by the lexico- 
graphically smallest set C{S) C S that generates S, i.e., cl(C(S')) = S. 
We call C{S) the canonical spanning set of the face S. The map C(-) is 
one-to-one; its inverse map is the closure map. 

The set C{S) can be computed efficiently as follows. For /c = 1,2 
set C{S) := S. For A; > 3, C{S) is computed iteratively: initialize C{S) 
with the set {si, S2}; at each iteration extend C{S) by the smallest Si such 
that d{C{S)) C c\{C{S) U {si}). Note that \C{S)\ < dim(5) + 1 < d + 1. 
Recall that we stored the vertex-facet incidences in the sorted sparse for- 
mat (see Section 2.2). Similarly to the method for computing closures, this 
computation can be performed in 0(a) steps, since just the intersections 
F({si}) n • • • n F({sj}), i = 1, . . . ,k, have to be computed iteratively. Then, 
C{S) is obtained as the set of those Sj for which the intersection becomes 
smaller. 
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We now explain the structure of the face tree. Its arcs are directed away 
from the root. They are labeled with vertex numbers, such that no two 
arcs leaving the same node have the same label and on every directed path 
in the tree the labels are increasing. Via the sets of labels on the paths 
from the root, the nodes of the tree correspond to the sorted sets C{S) for 
the faces S Q V that have been seen so far. In particular, the root node 
represents the face 0. Each node has a pointer to the corresponding node 
of C By construction, the depth of the tree is bounded by d + 1. 

Suppose we want to find the node £s corresponding to some face S CV 
in the part of C that we have already constructed or to assert that this 
face has not yet been seen. We first sort S (a subset of {!,..., n}) in- 
creasingly in 0(n) steps (by counting or bucket sort, see [2, Chap. 8]) and 
compute C{S) in 0(a) steps. Then, starting from the root, we proceed (as 
long as possible) downwards in the face tree along arcs labeled by the suc- 
cessive elements of C{S). Either we find an existing node in the tree which 
corresponds to S, or we have to introduce new labeled arcs (and nodes) into 
the tree until we have constructed a node representing S. 

In the latter case, it might be necessary to construct an entire new path 
in the tree. The definition of the canonical spanning sets C{S) ensures that 
all "intermediate nodes" on that path will correspond to canonical spanning 
sets of faces as well. Hence, the number of nodes in the face tree always 
will be bounded by (p, the total number of faces of the polytope. The faces 
represented by intermediate nodes will be seen later in the algorithm. Con- 
sequently, the corresponding pointers to C are set to nil for the meantime. 
Later in the algorithm, when we are searching for the face represented by 
such a tree-node for the first time, the nil-pointer will indicate that this face 
is not yet represented in C. The nil-pointer is then replaced by a pointer 
to a newly created node representing the face in C 

In any case, since the face tree has depth at most d+ 1 and the out-degree 
of each node is at most n, we need a total time of 0(n -|- a -|- (c? -|- 1) • n) = 
0(a) to either locate or create the tree-node representing a certain face. 

Lemma 3. Using the face tree, it is possible to locate or create the node 
in C representing a face in 0{a) steps (provided the vertex-facet incidence 
matrix is stored in the sorted sparse format). 

In the description given above, we have assumed that for each node in 
the face tree the out-arcs are stored in a list which is searched linearly for 
a certain label when walking down the tree. Of course, one can store the 
set of out-arcs in a balanced search tree (see, e.g., [2, Chap. 13]), allowing 



7 



to perform the search for a certain label in logarithmic time. After comput- 
ing C{S) for a face S (in 0(a) time), this allows to locate or create the node 
corresponding to S in the face tree in 0((d+ 1) • logn) steps. The total 
running time remains 0(a); nevertheless this might speed up the algorithm 
in practice. 

Instead of using the face tree, one can also store the faces in a bal- 
anced search tree. Again, the faces are represented by their canonical 
spanning sets, which are ordered lexicographically. Once C{S) is com- 
puted for a face S, searching S can be performed in 0{{d+ 1) • log(^) C 
0((d + 1) • min{n, m}) steps (since (p < 2™™{'^'™)). This yields the same 
total asymptotic running time, but searching the tree takes more (or the 
same) time compared to the variant of the face tree with balanced search 
trees at its nodes, since logn < min{ri, m}. 

2.5 The Analysis 

We summarize the algorithm in pseudo-code (Algorithm 1): 

Algorithm 1 Computing the face lattice of a polytope from its incidences 
1: Input: incidence matrix of a polytope P 
2: Output: Hasse diagram C of the face lattice of P 
3: initialize C and a face tree with £0 corresponding to the empty face 
4: initialize a set Q C {nodes of £} x {subsets of V} by {£0, 0) 
5: while Qy^ do 

6: choose some {£h, H) £ Q and remove it from Q 
7: compute the collection H ol all H{v), v \ H 
8: compute the set Q of minimal sets in H 
9: for each G eQ do 

10: locate/create the node £g corresponding to G in £ 
11: if -^G was newly created then 
12: add(£G,G)toQ 
13: end if 

14: add the arc {IrAg) to C 
15: end for 
16: end while 



Proposition 1. Algorithm 1 computes the Hasse diagram of the face lattice 
of a polytope P from its vertex-facet incidences in 0{n ■ a ■ (f) time. It can 
be implemented such that its space requirements (without output space) are 
bounded by 0{ip ■ n) . 
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Proof. Algorithm 1 works correctly by the discussion above. 

Step 7 can be performed in 0(n • a) steps by Lemma 1. Lemma 2 shows 
that we can execute Step 8 in O(n^) C 0(n • a) time. Hence, Steps 7 and 8 
in total contribute at most 0{n ■ a ■ (p) to the running time (since the while- 
loop is executed once per face). 

The body of the for-loop has to be executed for each of the 0{n ■ ip) arcs 
in the Hasse diagram C. Lemma 3 implies that each execution of the body 
of the for-loop can be performed in 0(a) steps. Thus, the claim on the 
running time follows. 

Since each node of the face tree corresponds to a face of P, the face tree 
has 0(99) nodes. Each label on an edge of the face tree needs at most 0(log n) 
bits, and we can estimate the space requirements of any of the (internal and 
external) pointers by 0(log(/?) C 0(min{n,m}). Hence, the face tree needs 
no more than 0{(p ■ min{n, m}) bits. 

The space required for the storage of Q is bounded by 0(^9 • n), if for 
each pair (iff, H) G Q the set H is stored as a hit set, i.e., the characteristic 
vector of C F is stored bit by bit. □ 

If m < n, then it is more efficient to apply Algorithm 1 to the incidences 
of the dual polytope, i.e., to the transpose of the incidence matrix. Of 
course, after the computations the roles of vertices and facets have to be 
exchanged again. This yields the main result of the paper. 

Theorem 1. The Hasse diagram of the face lattice of a polytope P can 
he computed from the vertex-facet incidences of P in 0(min{n, m} ■ a ■ p) 
time, where n is the numher of vertices, m is the numher of facets, a is 
the number of vertex-facet incidences, and (p is the total number of faces 
of P. The space requirements of the algorithm (without output space) can 
be bounded by 0{ip ■ min{n, m}). 

Whenever a new node representing a face G in the Hasse diagram C is 
constructed, we can label that node with the vertex set of G, the set of 
facets containing G, or with the dimension of G without (asymptotically) 
increasing the running time of the algorithm. The output, however, might 
become much larger with such labclings. For instance, labeling the Hasse 
diagram of the d-cube by vertex labels requires 17(4'^ • d) output storage 
space, while the Hasse diagram with facet labels needs only 0((i • S'^ • logd) 
space. 

In practice, the computation can be speeded up by exploiting that every 
vertex that is contained in a face G with H C. G and dim G = dim H -\- 1 
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must be contained in some facet which contains H. Thus, it suffices to 
consider only the sets H{v), v G (U/eF(if) ^({/})) \H m. Step 7. 

3 Extensions 

3.1 Simple or Simplicial Polytopes 

For a simple d-polytope P with n vertices, the above procedure can be 
implemented to run more efficiently. We have a = n - din this case. Prom 

the incidences (stored in the sorted sparse format), the graph G{P) of P 
(i.e., all one-dimensional faces) can be computed in time 0(n^ • d), since a 
pair of vertices forms an edge if and only if it is contained in d — \ common 
facets. 

After initialization with the vertices instead of (in Steps 3 and 4), 
Steps 7 and 8 can now be simplified. Consider an arbitrary vertex G H. 
For each neighbor v ^ H oi w in G{P), H{y) yields the other end node of 
an arc in the Hasse diagram; and each out-arc of H is produced this way. 
Thus, we can avoid constructing non-minimal faces in Step 7. Hence, Step 8 
can be skipped. The total running time for simple d-polytopes decreases 
to 0(d ■ a ■ Lp) (since the body of the for-loop is executed at most d ■ ip 
times) . 

The space complexity stays 0{ip ■ n) (see Proposition 1). It can, how- 
ever, be reduced to 0((/3 • m) (we have m < n for simple polytopes): in- 
stead of storing pairs {Ih-,H) in the set Q, we store the pairs {£h,F(H)), 
since |P(i?)| < m. Converting between H and F{H) can be performed 
in 0{a) steps and hence does not increase the asymptotic total running 
time. 

By duality, the same running times and space requirements can be 
achieved for simplicial polytopes. 

Similarly to the situation with general polytopes, for both simple and 
simplicial polytopes we can also output for each face its vertices, the facets 
containing it, or its dimension without (asymptotically) increasing the run- 
ning time. 

3.2 The fc-Skeleton 

A variant of Algorithm 1 computes the Hasse diagram of the fe-skeleton 
(all faces of dimension at most k) of a polytope P. One simply prevents 
the computation of faces of dimensions larger than k by not inserting any 
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{k — l)-face into the list Q. This leads to an 0(n • a ■ (p- ) time algorithm, 
where ip-'' is the number of faces of P of dimension at most k. 

3.3 Computing the "Hasse Diagram without Edges" 

If we only want to compute the faces of P together with their descriptions 
and dimensions (i.e., the "Hasse diagram without edges"), there exists a 
variant of Algorithm 1 with the same asymptotic running time, but signifi- 
cantly reduced space requirements. The difference is that no face tree is 
used, and the set Q is organized as a stack, i.e., the faces are investigated 
in a depth-first search manner. At each step, we take a face H from the 
stack, output it, and compute the set Q of {dimH + l)-faces containing H, 
like in Steps 7 and 8 of Algorithm 1. This needs time 0{n ■ a) for each H. 
The for-loop beginning at Step 10, including the search in the face tree, is 
replaced by an efficient way to decide which of the faces in Q is put onto 
the stack Q, such that every face appears exactly once on the stack during 
the algorithm. For this, we compute for each face G & Q a unique canonical 
facet H' of it. We put G onto the stack if and only ii H = H' . This ensures 
that each face is picked exactly once. 

We take H' as the closure of a set D{G), which is computed similar 
to the set C{G) of Section 2.4, except that we reject vertices which would 
produce G. More precisely, let G = {gi,g2, ■ ■ ■ ,gi}, with gi < §2 < ■ ■ ■ < 91- 
Initialize D(G) with and in each iteration extend D{G) by the smallest gi 
such that cl(L»(G)) C d{D{G) U {gi}) and cl{D{G) U {gi}) / G. After 
the computation, H' , the closure of D{G), clearly is a proper face of G. 
Moreover, it is a facet of G, since otherwise there exists a vertex g € G\H' , 
such that cl{H' U {g}) C G. But then g would have been included into D{G) 
when it was considered. Hence, D{G) is the lexicographically smallest subset 
of G which spans a facet of G. It can be computed in time 0(a), and hence, 
checking for all faces G G G whether H is the canonical facet D{G) of G can 
be performed in 0(ra • a) time. 

Altogether, this leads to an 0(n ■ a ■ p) time algorithm (see the proof of 
Proposition 1). The algorithm needs 0(n^ • d ■ logn) space for Q; since the 
depth of Q is at most d+l, there are never more than n- (d+l) sets on the 
stack, each of size at most n. Additionally, we need 0(a • log max{n, m}) 
space for storing the incidences in the sorted sparse format. Applying this 
method to the dual polytope, if necessary, we obtain an 0(min{n, m} ■ a ■ (p) 
time algorithm. 
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3.4 Oriented Matroids 



Algorithm 1 can be used for the enumeration of the elements of any atomic 
lattice provided a subroutine is available that computes the join of a set of 
atoms. For instance, this holds for every atomic and coatomic lattice if the 
atom-coatom incidences are given, because in this case one can compute the 
joins of atoms similarly to the case of face lattices of polytopes. 

In the following, we describe such an application of our algorithm to 
oriented matroids. The set of covectors of an oriented matroid with ground 
set {1, . . . , fc} is a subset of {— , 0, +}^ that satisfies certain axioms. We refer 
to Bjorner, Las Vergnas, Sturmfels, White, and Ziegler [1, Chap. 4] for the 
definitions and concepts that are relevant in the following. A specific, but 
illustrative, example arises from any finite set X of points in R,*^ as follows. 
For every linear functional ip G (1R°')* denote by SIGn((/?) G {—,0,+}'''- the 
vector whose component corresponding to x & X encodes the sign of (p{x). 
Then |sign((^) : ip G (R'^)*} is the set of covectors of an oriented mat- 
roid 0{X). 

For V, w G { — , 0, +}^ the separation set of V and w contains all indices i 
such that one of Vj, Wj is +, and the other one is — . The composition vow 
of V and W is defined by (vow)i := Vj if 7^ and (vow)j := Wj otherwise. 

We define a partial order ^ on {—,0,+}^, where V ^ W holds if and 
only if for all i we have Vj = or Vj = Wj. The ^-minimal elements among 
the nonzero covectors of an oriented matroid are called its cocircuits. If one 
adjoins an artificial maximal element 1 to the poset formed by the covectors 
of an oriented matroid (ordered by then one obtains its (big) face lattice. 

If, in the above example, X is the vertex set of a polytope P C R.*^, 
then the faces of P correspond to the positive covectors (i.e., the covectors 
with no component equal to — ) of 0{X). The facets of P correspond to the 
positive cocircuits of 0{X). The face lattice of P is anti-isomorphic to a 
sublattice of the face lattice of 0{X). 

The face lattice of an oriented matroid is atomic and coatomic; its atoms 
are the cocircuits, and its coatoms are called topes. Hence, we can compute 
its Hasse diagram from the abstract atom-coatom incidences as above. 

However, this is not the usual way to encode an oriented matroid. It is 
more common to specify an oriented matroid by its cocircuits. The join of 
two covectors simply is their composition, if their separation set is empty, 
or 1 otherwise. Such a composition can be computed in 0{k) steps, which 
enables us to compute the face lattice (efficiently) from its cocircuits by a 
variant of Algorithm 1. 

In Step 6, H now is a face of the oriented matroid, i.e., a covector. In 
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Step 7, one has to compute the joins of H with every cocircuit Y ^ H. Thus, 
Step 7 can be performed m. 0{n ■ k ■ (p) steps altogether (where ip is the total 
number of covectors and n is the number of cocircuits). We do not know 
any method to perform Step 8 faster than by pairwise comparisons, which 
take 0(n^ ■ k ■ (p) time in total. 

The face tree is organized similarly to the description in Section 2.4. 
One fixes an ordering Ci,...,C„ of the cocircuits. For a covector S let 
{ii, . . . , ir} {ii < ■ ■ ■ < ir) be the index set of cocircuits Ci- ^ S. Then we 
iteratively form the joins of Cj^, . . . , Cj^, and let C(S) consist of all those 
indices for which the "joins change." Computing C(S) from S takes 0{n ■ k) 
steps. Note that \C{S)\ < k. 

Using this modified face tree, a given covector S can now be searched 
in the same way as in the case of face lattices of polytopes. The depth of 
the face tree is bounded by k. Hence, location/creation of a covector can be 
performed in 0(n • k) time. The rest of the analysis is similar to the proof 
of Proposition 1. Thus, by this variant of Algorithm 1, the Hasse diagram 
of the face lattice of an oriented matroid can be computed in 0{n? ■ k ■ p) 
steps, requiring 0{ip ■ k) working space (since (p < 3^). 

Finschi [3] describes a different algorithm that computes the covectors of 
an oriented matroid from its cocircuits in 0(n • fc^ • (p) time. His algorithm, 
however, does not produce the edges of the Hasse diagram. 

The case where the topes (i.e., the :^-maximal covectors) of an oriented 
matroid are given is a bit different. Here, the number of faces is bounded 
by m^, where m is the number of topes. Hence, the size of the face lattice 
is polynomial in m. Pukuda, Saito, and Tamura [6] give an 0{k^ ■ m?) time 
algorithm for constructing the face lattice from the maximal covectors. 
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